Isotope shifts of the 4s 2 1 S'o — > 4s5p l P\ transition and 
hyperfine splitting of the 4s5p l P\ state in calcium 
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Using a technique based on production of ion Coulomb crystals, the isotope shifts of the 4s 2 
1 Sq — > 4s5p 1 Pi transition for all naturally occurring isotopes of calcium as well as the hyperfine 
splitting of the 4s5p 1 Pi state in 43 Ca have been measured. The field shift and specific mass shift 
coefficients as well as the hyperfine structure constants for 43 Ca have been derived from the data. 
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I. INTRODUCTION 

Isotope shifts and hyperfine splitting of optical transi- 
tions provide valuable information about atomic electron 
configurations and not at least properties of the nuclei. 
While from isotope shifts, nuclear charge distributions 
can be deduced (see, e.g., Refs 0, 0), from hyperfine 
splittings, nuclear spins, magnetic dipole moments, and 
electric quadrupole moments can be determined [3j. 

Although isotope shifts and hyperfine splittings al- 
ready have been measured for a large number of tran- 
sitions in neutral Ca (see Refs. |ji Lj l| and references 
therein) and in singly charged Ca + ions 0, IE the 
4s 2 1 So — * 4s5p 1 P\ transition considered in the present 
paper has hitherto not been measured. In contrast to 
previous experiments, the present measurements use ion 
Coulomb crystals as means to obtain the spectroscopic 
data. 



II. EXPERIMENTAL PROCEDURE 

In the experiments, we have used a frequency doubled 
cw laser system tuned to the 4s 2 1 Sq — » is5p 1 Pi tran- 
sition at a wavelength of A w 272 nm. The degree of 
excitation at a particular wavelength is determined from 
the ion production rate when excited neutral Ca atoms 
absorb yet another A ~ 272 nm photon directly from the 
4s5p 1 P\ state or from the metastable 4s3d 1 Dq state, 
populated through spontaneous emission as depicted in 
Fig. nja) [Tc| . As shown in the sketch of the experi- 
mental setup in Fig. Efa) , the spectroscopic laser beam 
at A w 272 nm crosses a well-collimated effusive ther- 
mal beam of calcium atoms at right angles in the center 
of a linear Paul trap ^jj. This geometry is chosen in 
order to obtain small Doppler shifts of the transition un- 
der study as well as to have a large capture efficiency 
of the ions produced. Since the atomic beam is derived 



from an oven containing natural metallic calcium, all the 
isotopes 40 Ca (96.9%), 42 Ca (0.647%), 43 Ca (0.135%), 
44 Ca (2.086%), 46 Ca (0.004%), and 48 Ca (0.187%) are 
present. The detection of the ions produced within a 
given exposure period is done by first trapping them in- 
dependently of the specific isotope, and then converting 
almost all ions into 40 Ca + through near-resonant elec- 
tron transfer collisions with atoms in the atomic beam 
containing 96.9% 40 Ca atoms. Finally, the number of 
trapped ions is counted with a near-100% efficiency by 
monitoring the fluorescence from the 40 Ca + ions when 
laser cooled into a Coulomb crystal 0,0]. Below, the 
key parts of the applied technique are discussed in detail. 
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FIG. 1: Simplified energy level diagrams of Ca (a) and Ca + 
(b) of relevance to the experiments. The straight arrows indi- 
cate laser driving transitions, while the wavy arrow symbol- 
izes spontaneous emission. In (a), the solid line beneath the 
hatched region indicates the ionization limit. 



A. Isotope and hyperfine selective ion production 

The laser system used for the spectroscopy consists of 
a tunable, single-mode ring CW dye laser operated at 544 
nm with an output power of a few hundreds mW which 
is frequency doubled to 272 nm by a /3-Barium Borate 
crystal placed in an external cavity. The output power 
at the desired 4s 2 1 So — > 4s5p 1 P\ transition wavelength 
is typically about ~ 10 mW, and the laser linewidth is 
below 1 MHz. This linewidth is much smaller than the 
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FIG. 2: (a) Schematics of the experimental setup with an inset 
showing a side view of the imaging system. For the optical el- 
ements in the figure, the following notation has been used: M, 
mirror; DM, dichroic mirror; PBS, polarization beam splitter; 
and A/2, half-wave-plate, (b) The linear Paul trap electrodes 
with applied rf voltages. The z axis defines the trap axis. 
The middle sections are 5.4 mm long and the end pieces are 
20 mm each. The diameter of the electrodes is 8.0 mm and 
the minimum distance to the trap axis is 3.5 mm. 



full width at half maximum(FWHM) Dopplcr width of 
~ 50 MHz due to the transverse velocity spread of the 
thermal atomic beam originating from an oven at a tem- 
perature of ~ 600°C. Hence, since the natural linewidth 
of the 4s 2 1 Sq — » 4s5p X P\ transition according to the 
NIST database [3 is known to be Y w 2 MHz (> 50% 
uncertainty) , the Doppler effect is the resolution limiting 
factor in the experiments. With isotope shifts typically 
of the order of GHz and hyperfine splittings of about ~ 
100 MHz of 43 Ca, not only can one selectively excite indi- 
vidual even isotopes (no nuclear spin) , but also the three 
hyperfine components of the 4s5p l P\ state in the case 
of 43 Ca. As mentioned above, ions are produced from 
excited calcium atoms either by absorbing yet another 
A s» 272 nm photon directly from the 4s5p x Pi state, or 
from the metastable 4s3d 1 Dq state populated through 
spontaneous decay. By exposing the atomic beam to the 
photoionization laser beam with a certain frequency for 
a fixed period of time (in our experiment controlled by a 
mechanical light shutter), the ion production will reflect 
the excitation degree of the 4s5p 1 P\ state of calcium 



atoms at this frequency. Counting the number of pro- 
duced ions as a function of laser frequency thus reveals 
the excitation spectrum. 



B. The ion detection scheme 

The central element of the ion detection is the linear 
Paul trap situated in a vacuum chamber which is op- 
erated at a pressure of ~ 10~ 10 Torr. As sketched in 
Fig- IHb), the trap consists of four sectioned cylindrical 
electrodes placed in a quadrupolc configuration. Con- 
finement of ions in the xy plane is obtained by an ef- 
fective radial harmonical potential created by applying 
a rf-potential of amplitude Urf and frequency £Irf to 
each of the electrodes in such a way that the voltage 
on the diagonally opposite electrode is in phase while 
the voltage of the neighbouring electrodes is 180° out 
of phase. Along the trap axis, the z-axis, confinement 
originates from a positive dc voltage Uec applied to all 
the eight end-pieces. With the trap dimensions used [see 
Fig- HM], for U RF ~ 200 V, ft = 2tt x 3.9 MHz and 
Uec ~ 5 V, the effective trap depth becomes ~ 1 eV. 
This ensures that ions produced in the trapping region 
are efficiently trapped despite the thermal velocities of 
the neutral atoms. 

When ions of a certain calcium isotope, e.g., 40 Ca + , 
are produced, they can be Doppler laser-cooled by driv- 
ing the 4s S1/2 — > 4p P1/2 and the 3d D 3 / 2 — ► 4p Pi/ 2 
transitions by laser light around 397 nm and 866 nm, 
respectively, as depicted in Fig. ^b). As the ions are 
cooled to sufficiently low temperatures (~ 10 mK), they 
arrange in a spatial ordered structure, often referred to 
as an ion Coulomb crystal. In our setup, these ion 
Coulomb crystals are imaged onto an image intensified 
Charge-coupled device camera by collecting fluorescence 
from the Doppler cooled ions. A typical picture of a 
Coulomb crystal containing ~ 10 000 40 Ca + ions is shown 
in Fig. |3| Here, one clearly observes the elliptical projec- 
tion of the crystal, which for our trapping potential is 
spheroidal with the trap axis (z axis) as the axis of rota- 
tional symmetry. Single-component Coulomb crystals in 
our trap have spatially uniform ion densities when viewed 
on length scales larger than the typical distance between 
neighboring ions (see Ref. ^H). The specific density is 
set by the trapping parameters. Hence, by determining 
the crystal volumes, e.g., by measuring the main axis of 
the projected ellipse, we can determine the total number 
of trapped ions. 

Generally, when the frequency of the spectroscopy laser 
is tuned between the resonances the various isotopes, two 
or more singly charged calcium isotope ions will be pro- 
duced and trapped. Such multispecies ion ensembles can 
be cooled into Coulomb crystals as well, with the lighter 
isotopes more tightly bound towards the trap axis than 
the heavier due to the dynamical confinement in the xy 
plane In Fig. a), an image of the fluorescence from 
laser cooled 40 Ca + ions in such a crystal is shown. The 



FIG. 3: Image of a pure Ca Coulomb crystal consisting of 
~ 10 000 ions. 
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crystal consists mainly of 40 Ca + - and 44 Ca + ions, where 
the nonfluorescing 44 Ca + ions are cooled only sympathet- 
ically through the Coulomb interaction with the directly 
laser cooled 40 Ca + ions [lit Il7j. The number of ions in 
the multicomponcnt crystal could principally be deter- 
mined by imaging the fluorescence from the individual 
isotopes. This is, however, technically very demanding, 
since each isotope ion would need its own set of laser fre- 
quencies due to isotope shifts of the transitions shown 
in Fig. mb). As an alternative route to quantify the 
total ion production, we choose to expose the multicom- 
ponent crystals similar to the one in Fig. 01a) to the 
thermal atomic beam of calcium until nearly all the ions 
(~ 96.9%) by near-resonant electron transfer collisions of 
the type 
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have been converted into 40 Ca + ions, and subsequently 
we measure the size of the now nearly pure 40 Ca + crys- 
tal. Since the energy difference between the two sides of 
Eq. QJ is just the isotope shift of the electronic levels, 
which is much smaller than the thermal energy of the 
atoms in the beam, no significant energy barrier exists 
for the process. Furthermore, since the kinetic energy 
of the produced 40 Ca + ions is more than an order of 
magnitude smaller than the trap potential depth, essen- 
tially 100% of the ions produced through Eq. stay 
trapped. In Fig. Eft)), the two-component crystal pre- 
sented in Fig. U a ) has been converted via the electron 
transfer process into a nearly pure Ca + crystal. 



III. SPECTROSCOPIC RESULTS 

A. The 4s5p 1 Pi excitation spectrum 

Figure [5] shows the relative ion production rate or, 
equivalently, the relative rate at which calcium atoms 
are being excited to the 1 P\ state, as a function of the 
detuning of the spectroscopy laser. These relative rates 



FIG. 4: Image of a 40 Ca + (fluorescing)- 44 Ca + (nonfluorescing) 
two-component crystal of ~ 1 500 ions before (a) and after 
(b) exposure to the atomic beam which leads to the electron 
transfer as described in the text. The white line indicates the 
outer boundary of the crystal. 



were obtained from knowledge of the volumes of the cre- 
ated Coulomb crystals, the laser exposure times, and 
the intensities of the spectroscopy laser. The spectrum 
contains resonances of all the naturally abundant iso- 
topes, including the one corresponding to 46 Ca which 
has a natural abundance of only 0.004%. The results 
presented in Fig. [S] have been obtained from three par- 
tially overlapping frequency scans. The three data se- 
ries were measured at slightly different oven tempera- 
tures: T oven = 612°C (A), T oven = 630°C (■), and 
T oven = 602°C (□). The relative rates of the A data 
and the I data have been normalized, so they share the 
same fitted maximum value and position at the 42 Ca res- 
onance peak. In the same way the data indicated by I 
and □ have been normalized to the three 43 Ca hyperfine 
peaks. 

The temperature was during each experiment kept sta- 
ble to better than ±2°C corresponding to a maximum 
uncertainty in the atomic density of about ±7% during 
the whole data accumulation time of ~ 40 min. Since the 
time to measure a single resonance is only a very small 
fraction of this time, any systematic errors in the reso- 
nance profiles due to fluctuating oven temperatures can 
be neglected. The Doppler broadening of the resonances, 
which only depends slightly on the atomic mass, is best 
found by fitting the 40 Ca resonance data to a Voigt pro- 
file. By doing this one obtains a total FWHM of ~ 54 ±5 
MHz dominated by the Doppler broadening. At the laser 
intensity which we apply (~ 100 mW/cm 2 ), we have as- 
sumed that neither the 4s 2 So — > As5p 1 P\ transition nor 
the following photoionization process is saturated. Con- 
sequently, the ion production rate is expected to be pro- 
portional to the laser intensity squared. We accounted 
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for fluctuations in the spectroscopic laser intensity dur- 
ing the data accumulation time by using this assumption. 
That this assumption is rather good is supported by not- 
ing that the relative measured resonance peak values in 
Fig.[S]are in good agreement with the natural abundances 
of the calcium isotopes, which are 40 Ca (96.9%), 42 Ca 
(0.647%), 43 Ca (0.135%), 44 Ca (2.086%), 46 Ca (0.004%), 
and 48 Ca (0.187%). 




Detuning relative to Ca resonance (GHz) 

FIG. 5: Semilogarithmic plot of the production rate of Ca + 
ions as a function of the frequency detuning of the spec- 
troscopy laser from the 40 Ca resonance. The data points 
presented are the results of three frequency scans made at 
slightly different temperatures [(a) T ovm = 612°C, (■) 
T oven = 630°C, and (□) T oven = 602° C]. The full line in- 
dicates a Voigt profile fit to the 40 Ca resonance within the 
frequency range ~ [-1300; 700] MHz. 



B. Isotope shifts 

Before deriving the isotope shifts from data like the 
ones presented in Fig. [S] the uncertainty of the individ- 
ual data points has to be estimated and taken into ac- 
count. The uncertainty is essentially due to three effects. 
First, due to intensity fluctuations of the laser intensity 
there will be a fluctuation in the number of ions pro- 
duced. Since for neighboring data points around a reso- 
nance this intensity fluctuation is maximally a few per- 
cent, the expected square dependence on the intensity of 
the ion production leads to an estimated uncertainty of 
~ 5%. Second, in determination of the Coulomb crystal 
volumes, both systematic and random uncertainties due 
to the measurements of the main axis of the elliptical 
projections of the crystals occur. Both these errors are 
of about ~ 5%, but since the systematic errors are equal 
for points symmetrically positioned around a resonance, 
it will have a small effect on the determination of the 
resonance frequencies. Hence, in our analysis, we have 
only accounted for the random errors of 5%. Finally, due 



to the finite number of ions produced within one mea- 
surement, we have included an uncertainty of the square 
root of the estimated number of ions. If all these un- 
certainties are accounted for in a weighted least-squares 
Gaussian fit, the uncertainty on the resonance position 
becomes lower than 2 MHz. 

A more critical error arises when we compare similar 
scans. Here we find that the measured resonance peak 
positions are associated with much larger uncertainties 
than the ones from the Gaussian fits. This additional 
error originates from local frequency drift of the spec- 
troscopy laser during a whole scan which typically lasts 
about 40 min. From a series of scans, we have found 
that this laser drift error leads to a rms uncertainty in 
the resonance frequencies of about 9 MHz. This last error 
source is the dominant one for the estimated errors of the 
isotope shifts given in Table |U In this table, the isotope 
shift for 43 Ca is found as the center of gravity of the three 
hyperfine components of 43 Ca which will be discussed in 
Sec. IIII CI In addition to these statistical errors, there 
is also a systematic uncertainty of ±1% arising from the 
calibration of the spectroscopy laser frequency scan to 
an optical spectrum analyzer with a known free-spectral 
range. 

For completeness, we have estimated that we might 
have introduced an unimportant error of a few hundred 
kHz due to the fact that we have determined the res- 
onance frequencies from single Gaussian fits instead as 
from a more realistic multipeak Voigt profile fit. 



Mass (A) Shift (MHz) 

42 967 ± 9 

43 1455 ± 9 

44 1879 ± 14 
46 2746 ± 16 
48 3528 ± 16 



TABLE I: Isotope shifts of the transition 4s 2 1 5" — > 4s5p 1 Pi 
in calcium derived from the experimental data like the ones 
presented in Fig. All shifts are with respect to the 40 Ca 
resonance. The errors stated represent one standard deviation 
originating from the statistical errors in the experiments. In 
addition, the shifts are subject to an overall linear scaling 
uncertainty of 1% due to our frequency scan calibration (see 
text). The isotope shift for 43 Ca is the center of gravity of 
the hyperfine components. 

The isotope shift for a given transition is usually de- 
scribed as a sum of the mass and the field shift in the 
following way Q 

6v AA '=M^^+FS(rY A ', (2) 

where M is the mass shift coefficient, A and A' denote 
the atomic masses of the two isotopes, F is the field shift 
coefficient, and S(r 2 ) AA is the difference in mean square 
nuclear charge radii between the isotopes. 
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The mass shift is usually written as a sum of the nor- 
mal mass shift (NMS) and the specific mass shift (SMS), 
which means that we can write the mass shift coefficient 
as M = M N ms + M S ms- Here the NMS coefficient 
is given by the simple expression M^ms = VQm e /m u , 
where vq is the transition frequency, m e is the electron 
mass, and m u is the atomic mass unit. The NMS origi- 
nates from the reduced mass correction for the electron, 
while the SMS comes from the change in the correlated 
motion of all the electrons (see, e.g., Ref. |18|L Subtrac- 
tion of the NMS from the total isotope-shift gives the 
residual isotope shift (RIS), 



8v 



AA' 
RIS 



Ms ms (A' - A) /A A' + F6(r ) 



2\AA' 



(3) 



Rewriting Eq. yj by multiplication by the factor 
AA'/(A' - A) leads to' 



A' A 
A'-A ( 



AA' 
'RIS 



= M, 



SMS 



F 



which shows that M$ms and F can be determined from 
a linear fit when the 5(r 2 ) AA 's are known. Using values 
of the 5(r 2 ) AA 's for calcium from Ref. 0, in Fig. we 
have plotted the left-hand side of Eq. Q using the values 
of Table CU as a function of [A' A/ (A' - A)]S(r 2 ) AA ' for 
the fixed value of A = 40. The SMS coefficient and field 
shift coefficient, obtained by weighted linear regression 
fit to the data points in Fig. are listed in Table ITU 
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FIG. 6: Plot of the values [A'40/(A' -40)]5vf/s derived from 



the measurements as function of [A' 40/ (A' - 40)]<5(r 2 ) 40A for 
the 4s 2 1 S'o — * 4s5p 1 Pi transition in Ca. The dashed line is 
a weighted linear regression fit to determine the coefficients 
Ms ms and F from the relation given in Eq. (J1J. 



Ms ms (GHz amu) 


F (MHz/fm 2 ) 


243 ± 3 ± 9 


-179 ±39 ±2 



TABLE II: The specific mass shift Msms and field shift F 
coefficients for the 4s 2 1 5'o — > 4s5p 4 Pi transitions of calcium 
derived from the the linear fit presented in Fig. E] The first 
stated uncertainty estimate is the one standard deviation ob- 
tained from the linear regression to the data of Fig. |5] The 
systematic error in the isotope shifts due to the frequency scan 
calibration is included as the second uncertainty estimate. 
For completeness the NMS coefficient is Mmms = 604.3 GHz 
amu. 



overlap with the nucleus compared to an s electron, so 
the main contribution to the field shifts is expected in 
both cases to come from the s electrons. 



C. 43 Ca hyperfine splitting 

The nuclear spin of I = 7/2 for 43 Ca leads to three 
hyperfine levels of the 1 P\ state with total spins F — 
5/2,7/2, and 9/2, respectively. The hyperfine structure 
(hfs) constants and isotope shift of 43 Ca are determined 
by fitting to the normal Casimir formula Q 



AE F =Av r 



A 



-C 



B%C{C + 1)-2I(I+1)J(J+1) 
1 (21 - 1)(2J - 1)1 J : 



(5) 



where C = F{F+ 1) -1(1+ 1) - J( J+ 1), Av cg is the iso- 
tope shift of the center of gravity of the hfs, and A and B 
are the magnetic dipole and electric quadrupole coupling 
constants, respectively. Several scans across the three 
hyperfine resonances were made to increase the level of 
confidence of the A and B constants. The hfs constants 
derived from these scans are summarized in Table IIIII 
while the center of gravity for 43 Ca has been given in 
Table UJ The small value of the B constant indicates 
that the magnetic dipole coupling has the most promi- 
nent contribution to lifting the degeneracy of the 4s5p 
1 P\ level for 43 Ca. Compared with the work of Ref. 0, 
where the hfs constants for the 4s4p l P% state in 43 Ca 
have been measured to be A = —15.54 ± 0.03 MHz and 
B = -3.48 ± 0.13 MHz, the 4s5p 1 P 1 state has an op- 
posite sign for the A constant and the B constant is the 
same order of magnitude or smaller. 



A (MHz) 


B (MHz) 


39.8 ±0.8 ±0.4 


-0.3±3±0.03 



The field shift coefficient for the 4s 2 l S 4s5p 
1 Pi transition ( Table HTjl is, within the stated error, al- 
most equal to the experimentally determined field shift 
coefficient for the 4s 2 1 Sq — > 4s4p 1 Pi transition of 
F = -175.8 ±1.2 MHz/fm 2 reported in Ref. @. This can 
be attributed to the fact that a p electron has negligible 



TABLE III: hfs constants for the 4s5p 1 Pi state in 43 Ca. The 
first stated uncertainties originate from statistical errors in de- 
termining the resonance positions, while the second account 
for the systematic errors due to frequency calibration uncer- 
tainty. 
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IV. CONCLUSION 

With a technique based on determining the number of 
ions produced, collected, and cooled into Coulomb crys- 
tals through resonant two-photon ionization, the isotope 
shifts of the 4s 2 1 Sq — > 4s5p 1 P\ transition for all nat- 
urally occurring isotopes of calcium as well as the hy- 
pcrfinc splitting of the 4s5p 1 P\ state in 43 Ca have been 
measured. The field shift and specific mass shift coef- 
ficients as well as the hyperfine structure constants for 
43 Ca have been derived from the data. Though in the 
presented experimental scheme, near-resonant charge ex- 
change collisions were used to gain information about the 
ion number, it is not a requirement for the technique for 



these to be used. Actually, the technique should be ap- 
plicable to many atomic species as long as some laser 
cooled ions are simultaneously trapped and can sympa- 
thetically cool the species of interest into a Coulomb crys- 
tal [lT], [3 \2m • From the spatial organization of the ob- 
servable laser cooled species, one can easily deduce the 
number of the, e.g., nonfluorescing ions of interest. 
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